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I discuss the recent implementation of medium-modified splitting functions within the HERWIG 
angular-ordered parton shower algorithm and present a few results on transverse momentum, 
energy and angular distributions. 
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Much work has been carried out through the years in order to theoretically describe jet quench- 
ing, one of the most striking observations of heavy-ion collisions, namely the suppression of par- 
ticle multiplicity at large transverse momentum (pr), as well as other related phenomena, such as 
the disappearance or distortion of the spectra of particles at opposite directions with respect to a 
reference one [p. A higher energy loss in a dense medium is a possible explanation of jet quench- 
ing: partons which in the vacuum are potentially produced at large pj can further radiate, which 
determines suppression of the high-pr spectrum and enhancement of the low-pr one. 

In order to describe the higher branching probability in a medium, a simple prescription con- 
sists in adding to the Altarelli-Parisi splitting functions a term depending on the medium properties 

P(z)^P(z)+AP(z,p 2 ,E,q,L). (1) 

In ([[J), p 2 is the branching-parton virtuality, E its energy, L the medium length, q the transport 
coefficient, defined as the average transverse momentum transferred from the medium to the parton 
per unity of free path. The correction AP(z,p 2 ,E,q,L) can be computed by means of the so-called 
BDMPS approximation [Q], namely multiple collinear radiation off static scattering centres, in the 
presence of a screened Coulomb potential. 

Such modified splitting functions can be implemented in Monte Carlo parton shower genera- 
tors, which can thus be used by the heavy-ion community to compare with data from RHIC and, 
ultimately, LHC experiments. Medium modifications have been included in the framework of the 
PYTHIA event generator [|5j], which orders cascades in virtuality, with an option to to veto branch- 
ings which do not satisfy angular ordering. This implementation is known as Q-PYTHIA [||, [7]]. 
More recently [|0, medium-modified splittings have also been implemented in the HERWIG show- 
ering algorithm [^], which, unlike PYTHIA, systematically includes angular ordering, thus leading 
to colour coherence in the large-Afc limit. Following herafter I shall present the main results 
yielded by medium-modified HERWIG parton showers. 

For the sake of comparison with the Q-PYTHIA results in [Q], I shall study the cascades 
initiated by a single gluon, whose energy will be fixed to E = 10 and 100 GeV, with the hadroniza- 
tion switched off. As discussed in [§], since this is not a standard HERWIG option, it was nec- 
essary to add a fictitious process and set by the hand the upper limit of the evolution variable 
Q 2 ~ E 2 (l — cos 6) to <2^ax = IE 2 , where E is the energy of the branching parton and 6 the emis- 
sion angle. Morover, care must be taken about the fact that the medium length varies along the 
shower: in fact, if Lq is the length for the first splitting, in a subsequent emission, a parton 'sees' 
an effective length L = Lq — 2zE /kj, z and kj being its energy fraction and transverse momentum, 
since it has travelled for a distance 2zE/k 2 , the so-called parton formation length, before radiating 
again. As L is not positive definite, whenever the medium length becomes negative, the shower 
will be vacuum-like. 

Hereafter, I shall consider the options of media with q = 1 and 10 GeV 2 /fm, Lq = 2 and 5 fm. 
For simplicity, such medium configurations will be labelled in terms of the so-called accumulated 
transverse momentum cjLq = 2, 5, 20 and 50 GeV 2 . 

A fundamental quantity in shower algorithms is the Sudakov form factor As(Q 2 ,Q 2 ), namely 
the probability of evolution between Q 2 and Q 2 with no intermediate branching. Fig. |I] presents the 
gluon Sudakov form factors As(Q 2 , <2max)> f° r E = 10 and 100 GeV, fixed L = Lq and the values 
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Figure 1: Gluon Sudakov form factors, in the vacuum (solid, black) and in media with qLo — 2 (dashes, 
blue), 5 (dots, red), 20 (dot-dashes, green) and 50 (magenta, solid) GeV 2 . The gluon energy is E = 10 GeV 
(a) and 100 GeV (b). 

Table 1: Average parton multiplicities in showers initiated by gluons of energy of 10 and 100 GeV, in the 
vacuum and in a medium with assigned values of qLo. 



E 


qL = 


qLo = 2 GeV 2 


qL Q = 5 GeV 2 


qL = 20 GeV 2 


qL Q = 50 GeV 2 


10 GeV 


2.56 


3.05 


4.14 


3.60 


4.56 


100 GeV 


6.95 


7.41 


8.79 


8.93 


11.70 



of q and Lq given above. From Fig. [I] one learns that medium-induced effects are quite large: 
the Sudakov form factors exhibit suppressions of a few orders of magnitude, with respect to the 
vacuum, corresponding to an enhancement of the branching probability, and thus of the radiative 
energy loss. 

As a further check of the higher emission probability, Table [I] quotes the the average parton 
multiplicity in the vacuum and in a dense medium. The enhancement factor runs from 20% (qLo = 
2 GeV 2 ) up to about 80% (c/L Q = 50 GeV 2 ) for E = 10 GeV, and between 7% and 70% for 100 GeV. 
In Figs. H-|] I present the differential multiplicity with respect to the parton transverse momentum 
(Pt), angle (6) and logathmic energy fraction, defined as ^ = \n(E/\p\), where E is the initiating- 
gluon energy and \p\ the modulus of the three-momentum of the partons in the cascade. All plotted 
distributions are normalized to unity. 

Overall, the medium-modified spectra present the features which one should expect and in 
agreement with the jet-quenching observations: suppression (enhancement) of the parton multi- 
plicity at large (small) pj, broader angular distributions, small-^ suppression. However, compared 
with the distributions presented in [^] and obtained with the Q-PYTHIA code, some differences 
can be noticed. As discussed in [§J, on average the emission probability in PYTHIA is larger than 
in HERWIG, as a consequence of the different evolution variables and the range wherein they are 
allowed to vary. For example, the vacuum spectra here obtained for E = 10 GeV exhibit a sharp 
peak at pj = 6 = E, = 0, corresponding to a fraction of events wherein the gluon evolves down 
to the infrared cutoff with no emission. This feature is less evident when running Q-PYTHIA [Q]. 
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Figure 2: Transverse momentum multiplicity in a shower initiated by a gluon of energy 10 GeV (a) and 
100 GeV (b), in the vacuum (solid, black) and in media with accumulated transverse momentum qLo = 2 
(dashes, blue), 5 (dots, red), 20 (dot-dashes, green) and 50 (solid, magenta) GeV 2 . 




As for the angular distributions, Fig. ||| shows that the multiplicity is negligible for 6 > 2; in fact, 
unlike PYTHIA, HERWIG presents a dead zone for large-angle radiation and we are studying the 
shower produced by a single parton, with no matrix-element matching. Ref. [||] also discussed 
medium-modified showers for fixed length L = Lq and reasonable results were obtained, namely 
stronger medium effects once L does not decrease throughout the cascade. 

In summary, I reviewed the main issues involved in the modification of the HERWIG angular- 
ordered shower algorithm for the sake of including the effects of a dense medium and possibly 
describing jet quenching. The release of such a code, which will be called Q-HERWIG |[T^], is 
currently in progress. Furthermore, in order to draw a final conclusion on the comparison with 
Q-PYTHIA, it will be compelling turning hadronization on and tune the two programs to the same 
data set. In this way, light will be shed even on the role played by the evolution variable in medium- 
modified cascades. 
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Figure 4: As in Figs. 2 and 3, but showing the logarithmic energy-fraction | . 
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